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Manipulation of sudden death of entanglement (ESD) of two qubits interacting with statistically 
uncorrelated thermal reservoirs is investigated. It is shown that for initially prepared X-states of 
the two qubits a simple (necessary and sufficient) criterion for ESD can be derived with the help 
of the Peres-Horodecki criterion. It is shown analytically that, in contrast to the zero-temperature 
case, at finite temperature of at least one of the reservoirs all initially prepared two-qubit X-states 
exhibit ESD. General conditions are derived under which ESD can be hastened, delayed, or averted. 
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I. INTRODUCTION 

Entanglement is a vital resource for quantum infor- 
mation processing Numerous processes relevant 
for quantum computation, quantum cryptography [2] 
or quantum teleportation pj] rely on entangled qubit 
states. Recently, Yu and Eberly [j, [E d, 0, found that 
reservoir-induced decay of single-qubit coherence can be 
slower than the corresponding decay of qubit entangle- 
ment. Both, abrupt and asymptotically gradual disap- 
pearance of entanglement were predicted in amplitude [5[ 
and phase damping channels Abrupt disappearance 
of entanglement in finite time was termed "entanglement 
sudden death" (ESD). Also recently, experimental evi- 
dence of ESD was reported for an optical setup @ and 
atomic ensembles Whereas first investigations on 
ESD concentrated on entangled two-qubit states, later it 
was also explored in a wider context and in higher dimen- 
sional Hilbert spaces [ul, E [S Q E d, d, [H • 

Clearly, ESD is a serious limiting factor for the use 
of entangled qubits in quantum information processing. 
It would seem important to stabilize quantum systems 
against this unwanted phenomenon. First studies in this 
context concentrated on changing initial states to more 
robust ones of the same degree of entanglement [20l . [2l| . 
Recently, we have addressed the practically relevant ques- 
tion whether it is possible to delay or even avert ESD by 
application of particularly chosen local unitary transfor- 
mations for a given initial state and a given open-system 
dynamics [22] . We demonstrated that this is indeed pos- 
sible for the special two-qubit system investigated first 
by Yu and Eberly B an d have found similar effects in 
qubit-qutrit systems |23j . 

In this paper, we generalize our previous results and 
investigate ESD of two qubits which are interacting with 
statistically independent (bosonic) reservoirs at finite 
temperatures. It is demonstrated that based on the 
Peres-Horodecki criterion [23, [2^ and on recent results 
of Huang and Zhu [1^ it is possible to develop system- 
atically a simple criterion capable of characterizing delay 
and avoidance of ESD of initially prepared two-qubit X- 
states in this open quantum system. With the help of 



this criterion it is proved that, in agreement with recent 
conjectures based on numerical case studies fl^l, all ini- 
tially prepared two-qubit X-states exhibit ESD if at least 
one of the statistically independent reservoirs is at finite 
temperature. However, if both reservoirs are at zero tem- 
perature, there are always some X-states for which ESD 
does not occur. Based on this criterion, we demonstrate 
that even at finite temperatures of the reservoirs it is pos- 
sible to hasten or to delay ESD. The characteristic time 
of ESD can be controlled by the time when appropriate 
local unitary operations are applied. However, unlike in 
the zero-temperature case, when at least one of the reser- 
voirs is at finite temperature it is not possible to avoid 
ESD completely by any choice of local unitary transfor- 
mations. 



II. OPEN-SYSTEM DYNAMICS OF TWO 
QUBITS COUPLED TO STATISTICALLY 
INDEPENDENT THERMAL RESERVOIRS 



In this section, we briefly summarize the basic equa- 
tions of motion governing our open quantum system of 
interest. In order to put the problem of ESD, its delay, 
and avoidance into perspective, let us consider two non- 
interacting qubits which are spatially well separated so 
that each of them interacts with its own thermal reser- 
voir. These two reservoirs are assumed to be statistically 
independent and are possibly also at different tempera- 
tures. The coupling of these two qubits to the reservoirs 
can originate physically from the coupling of two two- 
level atoms to the (resonant) modes of the electromag- 
netic radiation field, for example, with the local radiation 
field at thermal equilibrium. 

In the interaction picture and in the dipole- and rotat- 
ing wave approximation, the resulting equation of mo- 
tion of these two qubits is given by the master equation 
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FIG. 1: Schematic representation of the two non-interacting 
two-level atoms (qubits) A and B, initially prepared in an 
entangled state. Each of them interacts with its own local 
reservoir Ri and i?2. 
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For reservoirs representing the electromagnetic radiation 
field, TO and n denote the mean photon numbers of the 
local reservoirs coupling to qubits 1 and 2. The sponta- 
neous emission of atom i from its excited state to its 
ground state |0i) is described by the spontaneous decay 
rate 7^ and a\. are the corresponding raising {+) and low- 
ering (-) operators, i.e. (j\ — and cr'i — |Oi)(li|. 
The orthonormal atomic eigenstates |1) — \1) a ® 
|2) = |0)b, |3) = \Q)a ® |4) = \Q)a ® |0)b 
form the (computational) basis of the four dimensional 
Hilbert space of the two qubits. The derivation of Eq.([T]) 
also assumes the validity of the Born-Markov approxi- 
mation. The general solution valid for an arbitrary ini- 
tially prepared two-qubit state is given in Appendix A 
by Eq.jXT]). 

There is some numerical evidence [T3| that the pres- 
ence of non-zero mean thermal photon numbers in Eq. ([1]) 
may be responsible for ESD. However, systematic explo- 
ration of these phenomena which is capable of proving the 
sufficiency of non-zero photon numbers for ESD and of 
providing a systematic analytical understanding of ESD 
in zero-temperature reservoirs has been missing so far. 
Our main purpose is to close this gap. In particular, we 
shall develop a simple analytical criterion for ESD which 
allows a systematic understanding of ESD and its delay 
and avoidance for initially prepared entangled X-states. 



III. THE PERES-HORODECKI CRITERION 
AND ENTANGLEMENT SUDDEN DEATH 



In this section, we analyze conditions under which ini- 
tially prepared entangled two-qubit states evolving ac- 
cording to Eq.(IT]) exhibit ESD. Starting from the gen- 
eral time dependent solution of these equations ESD is 
analyzed with the help of the Peres-Horodecki criterion 
[2^ . [25! and with the help of all principal minors of the 
partially transposed time-dependent two-qubit quantum 
state. 

To check for separability of a quantum state, Peres [l^l 
and Horodeckis [23] developed a powerful necessary and 
sufficient condition valid for 2 (g) 2- and 2 (g) 3-systems. It 
states that in these systems a quantum state is separable 
if and only if its partially transposed density operator is 
also a valid quantum state. Based on this observation, a 
measure of entanglement, called negativity, was proposed 
[2^ . It is the sum of the absolute values of all the negative 
eigenvalues of the partially transposed density operator. 
For 2(8)2 systems, there can be at most one such negative 
eigenvalue [29|. 

Recently, Huang and Zhu [11] studied the Peres- 
Horodecki criterion by focussing on the principal mi- 
nors of the partially transposed density matrix. The 
principal minor [p^'^ {ijkl...)] of the partially transposed 
density operator p^^ is the determinant of the subma- 
trix A4{ijkl...) formed by the matrix elements of the 
i, j, k, I, ... rows and columns of the partial transpose p^^ , 
that is, 



/ fl^"^ fl^^ n^'^ 

PT „PT „PT „PT 



M{ijkl...) = 



PT 

Pkl 

PT 
PU . Plj 



Pkj ' 
PT 



Pjk ' Pjl 

PT PT 

Pkk ' Pjl ' 
PT PT 

Plk ' Pll ' 



(2) 



In general, if a matrix is positive semidefinite, then 
all its principal minors are non-negative, and vice versa 
[sot . Therefore, for an entangled two-qubit state p, the 
smallest principal minor of its partially transposed den- 
sity operator must be negative. For two-qubit states the 
non-negativity of the principal minors [p^'^ (2)], 

[pP^m, [p^^(4)], [pP^m, [pP^m, [p^^(24)], and 
[^■^-^(34)] is guaranteed already by the non- negativity of 
the original density matrix p. As a consequence, a general 
two-qubit state is entangled if and only if the minimum 
value of the remaining seven principal minors P{p^'^) is 
negative, that is. 



J 



Pip 



PT\ 



m^n{ [p^^ (14)], [p^^ (23)], [p^^ (123)], [p^^ (124)], [p^^ (134)], [p^^ (234)], [p^^ (1234)] } < 



(3) 
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Typically, the investigation of these seven principal mi- 
nors for general solutions p{t) of the density operator 
equation ([T]) is cumbersome. However, significant simpli- 
fications are possible for identical zero-temperature reser- 
voirs with m = n = 0, and 71 = 72 = 7 in Eq.([T]). Thus, 
Huang and Zhu [1^ could demonstrate that for asymp- 
totically long times with jt^ 1, the separability of 
p{t^) is determined by the initial state p. In particular, 
these authors showed that the matrix 



/ Pii P21 Pl3 

P12 pii + P22 pii 

P31 Pii Pll + P33 

P32 P31 + P42 P12 + P3i 




(4) 



which is defined in terms of the matrix elements of the 
initially prepared two-qubit quantum state p, determines 
the asymptotic separability of the two qubits and thus 
the presence or absence of ESD. 

If P{p) < 0, the two-qubit state p{t'^) is entangled and 
ESD does not occur. If P{p) > 0, the asymptotic quan- 
tum state p{t^) is separable and ESD takes place. In 
summary, Huang and Zhu .26] showed that for identical 
zero-temperature reservoirs the necessary and sufficient 
condition for ESD is given by 



Pip) > 0, 



(5) 



and is determined by the initially prepared two-qubit 
quantum state. 

For initially prepared X-states, that is, quantum states 
with P12 = pi3 = P24 — P3A = in the basis of Sec|TTl this 
condition can be simplified considerably. In Appendix B, 
all relevant seven principal minors of a given quantum 
state p are evaluated. From these expressions it is ap- 
parent that for initially prepared X-states, all the seven 
relevant principal minors of p of Eq. ^ are positive if and 
only if the two principal minors are positive, that is. 



[p(14)] , [p(23)] > 0. 



(6) 



Thus, in the case of identical zero-temperature reservoirs 
and for initially prepared X-states, only these two prin- 
cipal minors determine whether the time evolution of 
Eq.([T]) causes ESD or not. 

Even more general properties can be deduced from the 
explicit expressions for the principal minors of X-states as 
given in Appendix B. It is apparent from the general solu- 
tion of the density operator Eq. ([1]) as given in Appendix 
A that an initially prepared two-qubit X-state remains 
an X-state for all times. Combining this observation with 
the results of Appendix B leads to the general conclusion 
that initially prepared two-qubit X-states exhibit ESD if 
and only if at asymptotic times the principal minors 
are both positive, i.e. 



[p(i4)](t-) , wmt'oo) > 0. 



(7) 



This general result holds even if the two qubits interact 
with statistically independent reservoirs at finite temper- 
atures so that their dynamics are described by Eq.([T|). 



IV. TWO-QUBIT X-STATES AND QUANTUM 
CONTROL OF ESD 



In this section, we specialize our discussion of ESD 
to initially prepared arbitrary two-qubit X-states. De- 
lay and avoidance of ESD of initially prepared two- 
qubit X-states coupled to statistically independent zero- 
temperature reservoirs is discussed in subsection A. In 
subsection B these results are generalized to reservoirs at 
finite temperatures. In particular, it is proved that if at 
least one of the reservoirs has nonzero temperature, all 
initially prepared X-states exhibit ESD. 

Let us first of all briefly summarize some basic prop- 
erties of X-states. The density matrix of a two-qubit 
X-state is of the general form 



Px 



^ Pll pi4 

P22 P23 

P32 P33 

V P41 P44 



(8) 



P12 



Pi3 = P24 = P34 = 0. In particular. 



Werner states [31| are special cases of such X-states 
and some aspects of their ESD have already been dis- 
cussed 0, HI m [Hi. Eq. © describes a quantum 
state provided the unit trace and positivity conditions 

YA=iPti = li P22P33 > |p23p, and P11P44 > |pi4p are 
fulfilled. X-states are entangled if and only if either 
P22P33 < |pi4p or /011P44 < |p23p- Both conditions can- 
not hold simultaneously [29^. 



A. Delaying and avoiding ESD in statistically 
independent zero-temperature reservoirs 

As discussed in the previous section, for zero- 
temperature reservoirs the criterion for ESD is given by 
Eq.® which together with the results of Appendix B 
yields the necessary and sufficient conditions for ESD: 



[p(14)] = Pll - |P23P > 0, 

1/5(23)] = (pii+p22)(pii 



P33j 



|P14| 



> 0. 

(9) 



Depending on the degree of entanglement of the initially 
prepared two-qubit state, two different cases can be dis- 
tinguished. 

Case 1: For initially prepared entangled two-qubit 
states fulfiUing the condition. 



P11P44 < IP23I 



(10) 



the analytical expression for the negativity of the quan- 
tum state p(t) satisfying Eq.(IT|) for 71 = 72 = 7 and 
TO = n = is given by 
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Ni{px{t))=max[Q, ^ Fip, pu)^ - 4p^{piiFip, pu) - pi, - IP23P) - Fip,pu)] 



(11) 



with F{p, pu) = (1 - 2p + 2p^)pn + (1 - P){P22 + P33) + 
P44, p — exp(— 7t). For any initially entangled two-qubit 
state p, Eq. pUj) implies P22P33 > |pi4p so that one of 
the conditions of Eq. ([9]) is satisfied. Thus, provided also 
the other condition, namely pu > |p23pj is satisfied, an 
initially prepared entangled two-qubit state exhibits BSD 
and its negativity becomes zero at a finite time, say ti. 

Provided both conditions of Eq. ^ are fulfilled, ESD 
can be delayed or even avoided by local unitary opera- 
tions acting on the two qubits involved. In particular, 
let us concentrate on local unitary operations which ex- 
change the density matrix elements pii(i) and p^it) of 
the quantum state at a time t < ti in such a way that 
their product, i.e. pii{t + 0)p44{t + 0), remains constant 
but that the condition p44{t + 0) > |p23(i + 0)P is vio- 
lated. According to Eq.®, in such a case ESD will be 
avoided. p44{t) is the probability of finding both qubits 
in their ground states. Thus, as a consequence of the 
dynamics of Eq.(IT]), the density matrix element p44(t) 
increases monotonically. There will be a limiting time 
tjio for any possible switching of these matrix elements 
for which ESD can still be avoided. If the local operation 
is applied after this limiting time, ESD may possibly be 
delayed but it is unavoidable. 

This simple consequence of the criterion of Eq . (I9|) ex- 
plains recent numerical work on this problem [22i] . In 
fact, operations of this type can avoid ESD for any ini- 
tially prepared two-qubit X-state provided they are ap- 
plied at a time t < tsw In particular, this applies to the 
subset of Werner states with pi4 = which are mixtures 
of a singlet state with probability a and a completely un- 
polarized (chaotic) state. These Werner states exhibit 
ESD in the parameter range a e [1/3, (—1 + \/5)/2) 
(20I I where pii{t) > |/923(i)P while entanglement de- 
cays asymptotically for values of a in the range a £ 
((— 1 + \/5)/2, 1] which corresponds to pii{t) < \p23{t)\^- 

Let us now deal with the question of which unitary 
transformations can achieve such a switch between P44{t) 
and pii(t). The most general 2x2 unitary matrix acting 
on a qubit is given by 



ditions 



U{2) 



sin(6l)e^(^+'^) 



sin(6')e*("-'^) 



(12) 



which is a linear superposition of the Pauli matrices a^, 
ay, az, and the Identity matrix (Tq- Exchanging the ma- 
trix elements pu and P44 can be achieved by applying two 
appropriately chosen unitaries C/4 and Ub of the form of 
Eg. p^ on qubits A and B a.t a suitably chosen time, say 
t, so that the X-state px{t) is transformed into another 
X-state {Ua 8) UB)px{t){u\ (g) C/jj), for example. 

The most general local unitary operations transform- 
ing an arbitrary X-state into another one fulfill the con- 



sin(26lA) = sin(26'B) =0 — > 
9a = rATT/2, 9b = rB7T/2, 



(13) 



with rA,rB G X-state preserving local unitary trans- 
formations with even values of ta and do not have any 
significant effect on the density matrix elements except 
multiplying puit) by a constant phase factor. Odd val- 
ues of VA and vb serve the purpose of exchanging pii(t) 
and P44{t). For any odd value of = r^, for example, 
the corresponding unitary two-qubit operator is given by 



U 







g2i(/34 























(14) 



A case in which such a X-state-preserving local unitary 
transformation is applied only onto qubit B can be de- 
scribed by parameters 9a = cha = 13a = = 0, for ex- 
ample. They lead to the transformations pii{t) P22{t), 
Psait) P44(i)i and Pi4{t) P23{t)- In view of the char- 
acteristic time evolution of p22{t) in zero-temperature 
reservoirs (compare with Appendix A) and the criterion 
of Eq.® this implies that such a switch of matrix ele- 
ments may delay ESD but it cannot be avoided. 

Case 2: For initially prepared entangled two-qubit X- 
states satisfying the alternative condition 



/522P33 < |Pl4| 



(15) 



the analytical expression for the negativity of the result- 
ing quantum state px (t) is given by 



N2[px{t))= max [0, pW{p22 - Pszf + 4|/0i4p 

-(p22+P33)-(2-2p)pn)]. (16) 

Eq. (fT5| implies piip44 > |p23p- Thus, the first condition 
of Eq. ^ is always satisfied so that ESD occurs whenever 
also the second condition is satisfied. The simplest case 
arises for P23 = where the initially prepared state is 
a Werner state, i.e. an incoherent mixture of a triplet 
state with probability a and the completely unpolarized 
state. ESD takes place in the parameter regime a € 
[1/3,1) where both conditions are satisfied during the 
time evolution. In the case of an initially prepared Bell 
state, i.e. for a — 1, the second condition of Eq. ^ fails 
and entanglement decays asymptotically. 

As discussed above, the first condition of Eq. ([9]) is al- 
ways fulfilled in the cases considered here so that ESD 
takes place always except in the particular case of an 
initially prepared Bell state which fulfills the condition 
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(pii(t) + P22W)(PiiW + P33(<)) = IpmWP- As a con- 
sequence any switch capable of exchanging P2z{t) and 
Pi4(t) will be sufficient to avoid or delay ESD. Such a 
switch can be implemented by a local unitary X-state- 
preserving transformation acting on qubit A or B only. 
As a result the second condition of Eq. ([9]) remains always 
true, while the validity of the first condition depends on 
the choice of the switching time t. If P33(t) > |/9i4(t)P 
ESD is unavoidable. However, for all switching times vi- 
olating this condition ESD is averted completely. As an 
example, let us consider Werner states with P23(0) — 0. 
In the parameter range a G [1/3, (-l + \/5)/2) the condi- 
tion Pz3,{t) > \pu{t)\'^ is fulfilled for all times so that ESD 
takes place. In the parameter regime a G ((— 1-|-\/5)/2, 1] 
this condition is violated so that entanglement decays 
asymptotically. 

B. Controlling ESD in statistically independent 
finite-temperature reservoirs 

According to Eq.®, ESD takes place if and only if the 
two principal minors 

[p^^(14)] = m\m + If + e-(2"+i)7*[j;^(i)] 

_^g-2(2m+l)7tj^(2)j g-3(2m+l)7tj^(3)j 
_^p-4(2m+l)7tj^(4)j^ 

[p^^(23)] = m\m + if + e-(2™+i)7t[G(i)] 

_^g-2(2m+l)7tjg(2)j ^ g-3(2m-Hl)7tj^(3)j 

_^g-4(2m+l)7t[G(4)], (17) 

are positive in the limit of very long interaction times. 
For simplicity, we have taken m = n and 71 = 72 = 7 
in Eq. pT]) . The quantities F^*) and G^'^ are functions 
of m and of the initial matrix elements of the initially 
prepared quantum state. Their explicit forms are given 
in Appendix C. Furthermore, the general expressions for 
[/0'^"^(14)] and [p'^"^(23)] are provided in Appendix B. 
For sufficiently long times, say i > t~ and for m > 0, 
factors of the form e~'2'"+^'''''*°° are exponentially small 
and therefore both [p^-^(14)] and [p'^"^(23)] are positive. 
Analogously, one can show that for unequal values of the 
mean photon numbers m and n, both minors are positive 
if and only if at least one of these mean photon numbers 
is not equal to zero. Hence, we arrive at the central re- 
sult that if one of the (photon) reservoirs is at nonzero 
temperature all initially prepared X-states exhibit ESD. 

As ESD is unavoidable in these cases, it may be use- 
ful at least to delay it. Indeed this can be achieved for 
all possible X-states [1, [l^l • For the sake of demonstra- 
tion let us consider the particular example of an initially 
prepared entangled state of the form 

P=^(|l,l>(l,l| + 2|*)(vl'|) (18) 



with I*) = (|0, 1) ± |1, 0))/\/2. It is known that this en- 
tangled state, while interacting with a vacuum reservoir, 
looses its entanglement at « 0.5348 [13]. However, 
while interacting with reservoirs at finite temperatures, 
the time of sudden death for this initial state depends 
on the values of m and n. The solution of Eq. ([T]) for 
the input state of Eq. (fT8|) can be obtained easily us- 
ing the general solution given by Eq. (jAip . After taking 
the partial transpose of the resulting density operator 
it is possible to obtain an analytical expression for the 
negativity of the quantum state at any time t. Setting 
m = n = 0.1, for example, we observe that ESD oc- 
curs at time Iesd ^ 0.41 15/7. Depending on the time 
when local unitary transformations are applied to qubits 
A and B, ESD can be speeded up or delayed for some 
finite time. 
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FIG. 2: Dependence of the time t^nd of ESD on the switching 
time tsw for m = n = 0.1: The X-state-preserving local uni- 
tary transformations switch the density matrix elements pn 
and P44 in Eq. (|18p . Starting on the right at switching time 
tA ~ 0.4115/7 this dependence exhibits a broad and small dip 
before rising to the maximum possible time tend ~ 0.9817/7. 

Fig. ^ displays the time tend at which ESD takes 
place and its dependence on the time of switching 
The earlier appropriate local unitary transformations are 
applied, the more ESD is delayed. However, typically 
such a delay is possible only for a certain range of switch- 
ing times tsw, such as tgu, < ts = 0.279/7 in Fig. ([2]). 
Eventually ESD is unavoidable. In the case considered 
in Fig. (21), it takes place at t^nd ~ 0.9817/7. 

In Fig. ([3]) the relation between tend and tg^ is depicted 
for mean thermal photon numbers m = n ^ 0.01. In this 
case ESD occurs at tESD ~ 0.5172/7. If the switch is 
applied before tA ~ 0.5172/7, ESD is hastened. Any 
switch made before ts ~ 0.2877/7 delays ESD up to the 
maximum possible time tend ~ 2.7087/7. This larger 
delay in comparison with the case considered in Fig. ((21) 
is due to smaller values of the mean photon numbers. 

Two recent papers that have appeared since comple- 
tion of our work examine related topics. One studies two 
harmonic oscillators coupled to a common environment 
also modelled as oscillators [34]. ESD and entanglement 
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FIG. 3: Dependence of the time tend of ESD on the switching 
time tsw for m = n = 0.01: The local unitary transformation 
are the same as in Fig. ([2}. For switching times below Ia ~ 
0.5172/7 this dependence has a broad and small dip before 
rising to the maximum possible time tend ~ 2.7087/7. 



basis introduced in Sec. II. These solutions are given by 



to , n/o rrrimn +m[{n + l)pn+ 
[Zm + l)[Zn + Ij 

-nip22 - P33 + P44 ) ]e-(2"+i)72t + n[{m+ 1) 

pii + (m + l)p22 - m(p33 + Pii) ]e-(2™+i)7it 

+ [ (to + l)(n + l)/9ll - TOP33 - n{p22 + TO/922 
+TOP33 - mp44 ) ]e-[(2™+l)7i + (2n+l)72]t}^ 



^ (2TO+l)(2n+l) ^"^^" + ^^'^"[("+^^^" 

+P33 - n(p22 - P33 + P44)]e-(2"+l)^^* + 

(n + 1) [ (1 + to) pii + (to + 1) P22 - w 
x(p33 + P44)]e-^'™+'^^^* + h(m+l) 
X (n + 1) pii + to P33 + n ( (to + 1) P22 + 

TOP33 - mp44 ) ] e-[(2™+l)7i + (2n+l)72]t}^ 



revival are examined as well as non-Markovian environ- 
ments. The other studies two interacting qubits in a 
magnetic field and a thermal Markovian environment, 
presenting nonmonotonic relaxation rates as functions of 
the magnetic field and temperature (35l |. 



V. CONCLUSIONS 

A criterion has been presented characterizing the con- 
ditions which lead to ESD of X-states of two qubits cou- 
pled to statistically independent reservoirs at finite tem- 
peratures. Based on this criterion, we have presented an 
analytical description of ESD of X-states and its delaying 
or its avoidance by local unitary actions. We have proved 
that if at least one of the reservoirs is at finite tempera- 
ture, all X-states exhibit ESD. Thus, in these cases ESD 
can only be delayed but not averted. Preliminary studies 
[2^ indicate that similar results also hold for qubit-qutrit 
systems. 
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-{n{m + 1) + (m + 1) 



(2to + l)(2n+ 1) 

X [ (n + 1) pii + P33 - n(p22 - P33 + P44 ) ] 

xg-(2„+i)72t ^n[im + l)pn + {m+ 1) P22 
-TO (p33 + P44 ) ] e-(2"+i)''i* + [-(to + 1) 
X (71 + 1) pii + TO P33 -I- n ( (to + 1) P22 

+TO P33 - m P44 ) ] e-[(2™+l)7i + (2n+l)7.]t}^ 



Pi2{t) - :^^{m(pi2 + P34)e^('"+''^^* + [(m-f 1) 
Zm + 1 



xpi2 - TOP34J e 2 



[2(2m+l)7i+(2n+l)72]t 



Prsit) = ^^{n{pu + p2i)e^^'"'+'^'''' + [{n+l) 

XP13 - np24] eT^[(2™+l)7i+2(2n+l)72]t |^ 
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APPENDIX A 



P24(t) = ^;^{(n+l)(pi3 + P24)e^("^'^^* + 

[np24 - (n + l)pi3] e^[(2™+l)7i+2(2«+l)72]t }^ 



P34(i) = ^^{(m + l)(pi2+P34)e^(2"+i)^^* + 

[TOP34 - (m + 1)P12] eT^[2(2™+l)7i + (2n+l)72]t 



Solutions of Eq.([T|) for some special initially prepared , , _[(„i+i)..yj+(„+i)^2]t 

X-states have been provided in the Appendix of Ref. ■ ~ ^^"^ ^ ' 

In this Appendix we provide the most general solutions of 

Eq. ([T]) for any initially prepared two-qubit quantum state ^ ^ +( + ^ ) ]t 

p with density matrix elements py in the computational P23{t) = P23 e ^ 2 72 ^ (Al) 
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APPENDIX B 



For X-states, i.e. quantum states with pi2 — pis = 
P24 = P34 = Oi the dependence of aU seven principal- 
minors of Eq. ^ on [^■''-^(14)] and on [p-^-^(23)] is given 

by 



[p^^(23)] = 



Pii P23 

P32 P44 



P22 Pl4 
P41 P33 



Pll PiA 



iP23r 



P22 P33 - |Pl4p 



[p^^(123)] = 


Pll 





P22 




Pl4 


= p^^ [p^^(23)] , 







P41 


P33 




[p^^(124)] = 


Pll 





P22 


P23 



= P22 [P^^(14)] , 




P32 





P44 




[p^^(134)] = 


Pll 






P33 


P23 




= P33 [P^^(14)] , 




P32 





P44 




[p^^(234)] = 


P22 
P41 




Pl4 
P33 






P44 


= P44 [P^^(23)] , 



(1234)] 



Pll P23 

P22 Pl4 

P41 P33 

P32 P44 

[p^^(14)] [p^^(23)] 



= -2pJi (m + 1)^ + (to + l)pii [(2m2 + m- 1) 
X (p22 + P33) + 2 mp44 ] + m [ (to + 1) pl^ 
+ {2{m, + 1) P33 - TO (2to + 3) pii) P22 + 
(to + 1) + 2 - to (2to + 3) P33 P44 ], 



F^"') = [m^ (pii - P22 - P33 +P44) + 

TO (2 Pll - P22 - P33) + Pii]^- (CI) 

Similarly, the expressions for G*^*-* in Eq. (|17p are given 



by 



= TO (m + 1) [ (2to + 1) Pll + P22 + P33 — 2 TO P44 ] , 



-2to^ [2p22 - 4p33 P22 - P22 + 2p33 + Pll - P33 
-|-8|pi4p + P44] - 2to^ [4p22 - 8p33 P22 - 2p22 

+4^33 + 3pii - 2p33 + 16|pi4p + P44 ] + rn^ 

[pll - 2(P44 + 3) Pll - 6p^2 - 6P33 + P44 + 2P22 
+ 12 P22 P33 + 2 P33 - 24 |pi4p ] + TO [2 + (p22 
+P33 - 2 P44 - 2) Pll - 2 P22 - 2 P33 + 4 P22 P33 
-8|pi4p - P22 P44 - P33 P44 ] + Pll + P22 P33 
+P1I (P22 + P33) - |Pl4p, 



(Bl) 



APPENDIX C 

The expressions for i^^'^ in Eq. p?)) are given by 

= to(to + 1)[(2to+1)pii+P22+P33-2top44], 



= -2TO4[2p24~p44 + p22 + 8|p23p + P33] + 

2to^ [-2p44 + 2p33P44 + P44 - 16 |p23p - 2 P33 
+ 2 P22 (P44 - 1) ] - TO^ [ P22 + (2 P33 - 4 P44 
+ 3)p22 + P33 + 24 |p23|^ + 3p33 - 4p33 P44 
-P44 ] - 4 TO (to + 1)^ Pll - TO [p22 + 2 P33 P22 
+P22 + P33 + 8 |P23P + P33 ] - |P23P + (m + if 
Pll[(8p44 + 2) TO^ + (-4p22 - 4p33 + 2)to + 1 ], 



G(3) 



-2 pI^ (to, + \f + (to, + 1) Pll [ (2 TO^ + TO - 1) 

X (p22 + P33) + 2 TO P44 ] + m [ (to + 1) P22 + ( 2 

(to + 1) P33 - TO (2to + 3) P44 ) P22 + (to + 1) P33 
+ 2 TO^ P44 - TO (2to + 3) P33 P44 ] , 



G^"^^ = [m.^(pii - P22 - P33 + P44) + 

TO(2pii - P22 - P33) + pii]^. (C2) 
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